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Abstract 

We have studied responses to applied external forces of the quantum (N$ + Nb) model for 
-Afe-body interacting harmonic oscillator (HO) system subjected to iV^-body HO bath, by using 
canonical transformations combined with Husimi's method for a driven quantum HO [K. Husimi, 
Prog. Theor. Phys. 9, 381 (1953)]. It has been shown that the response to a uniform force 
expressed by the Hamiltonian: Hf = —f(t) J2k=i Q k ^ s generally not proportional to N$ except 
for no system-bath couplings, where f(t) expresses its time dependence and Qk denotes a position 
operator of fcth particle of the system. We have calculated also the response to a space- and time- 
dependent force expressed by Hf = —f(t)Y2k=iQk e l2nku ^ Ns , where the wavevector u is u = 
and u = —Ng/2 for uniform and staggered forces, respectively. The spatial correlation T m for a 
pair of positions of Qk and Qk+m has been studied as functions of N$ and the temperature. Our 
calculations have indicated an importance of taking account of finite N$ in studying quantum open 
systems which generally include arbitrary numbers of particles. 
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I. INTRODUCTION 



In recent years, there has been considerable interest in open small systems, whose physical 



properties have been studied both by experimental and theoretical methods We may 
prepare desired small systems by advanced new techniques. Theoretical studies of open 
systems have been made with the use of the Caldeira-Leggett (CL) type models 2j-|6|. 
CL-type models have been extensively studied by using various methods such as quantum 
Langevin equation and master equation Q]. The original CL model considers a system of 
a single particle (N$ = 1) which is subjected to a bath consisting of infinite numbers of 
uncoupled harmonic oscillators (HOs) (N B = oo). Recent studies with the CL model have 
tried to go beyond this restriction on N$ and N B . References {isl have employed the CL 
model with N$ = 1 and N B — 1 — 800 for studies of properties of small system cou pled 
to finite bath. CL-type models with N$ = 2 and N B = oo have been investigated jiol. 



Reference [12J discusses the master equation of arbitrary Ns system coupled to an arbitrary 
N B bath. In our previous study [3], we have adopted the (N s + N B ) model for iV^-body 
system subjected to iV^-body bath in order to calculate energy distributions of a system, 
which show intrigue properties as functions of N s , N B and a system-bath coupling. 

In adopting the CL-type model, we have implicitly assumed that physical quantities such 
as the energy and specific heat of a system with finite Ns (> 1) are given as N$ times of 
results of a system with Ns = 1. Our recent calculation [14(, however, has pointed out that 
it is generally not the case because the system specific heat, Cs(T; Ns, N B ), of the (Ns + N B ) 
model at temperature T is given by 



C S {T- N s , N B ) ^ N s C S {T; 1, N B ), (1) 

except for no system-bath couplings and/or in the high-temperature limit. Furthermore it 
has been shown that the low-temperature specific heat may be negative for finite Ns with 
a strong system-bath coupling [Ljj]. This is in contrast with Refs. 15, 16, 18] showing a 
non-negative system specific heat for HO system in CL-type models with (Ns, N B ) = (1, 1) 
and (l,oo). These results imply that we should explicitly take into account finite Ns in 
studying open systems which may generally include arbitrary numbers of particles. It is 
interesting and necessary to study responses to applied external forces of the (Ns + N B ) 
model, which is the purpose of the present paper. Responses of the CL models have been 
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mostly made for infinite baths for which Ohmic and Drude models are adopted (e.g., Ref. 
191 ] ) [6j. In this study, we employ the identical- frequency model for finite baths 14]. 

The paper is organized as follows. In Sec. II, we briefly explain the (JVg + Nb) model 



13l . Il4j , to which we apply the canonical transformations in order to obtain the diagonalized 



iamiltonian including external forces. By using Husimi's method for a driven quantum HO 



2l| . we calculate the response of the open HO system to sinusoidal and step forces. In 



Sec. Ill, we calculate also the response to space- and time-dependent forces. The spatial 
correlation T m between positions of two particles separated by a distance m is evaluated. 
The final Sec. IV is devoted to our conclusion. 



II. THE (N s + N B ) MODEL 



A. Quantum Langevin equation 



We consider the (N$ + Nb) model in which the a one-dimensional iV^-body system (Hs) is 
subjected to an iV^-body bath (Hb) by the interaction (Hi) j3,[l4]. The total Hamiltonian 



is assumed to be given by 



H = H s + H B + Hi, 



with 
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k=l n=l 

Here Pk (p n ) and Qk (q n ) express the momentum and position operators, respectively, of a 
HO with a mass of M (m) in the system (bath), D and K denote force constants in the 
system, uj n is the oscillator frequency of the bath, Ck n is a system-bath coupling and fit) 
stands for an applied force. Operators satisfy commutation relations, 



[Qk, Pi. 



ih5, 



kC- 



[q n ,p ri 



ih5 n 



[Qk, Qi) = [Pk, Pi] = [q n , q m ) = [p n ,p m ] = o. (7) 



Equation ([3]) expresses the interacting HO system for D^O and K ^ 0. In the limiting case 
of K = 0, the system consists of a collection of uncoupled (independent) HOs. The system 
is subjected to a bath consisting of a collection of uncoupled HOs with oscillator frequencies 
of {u n }. 

In conventional approaches to the quantum system-plus-bath model, we obtain equations 
of motion for Q k and q n , employing the Heisenberg equation, 

ihO = [0,H], (8) 

where O expresses an arbitrary operator and a dot stands for a derivative with respect of 
time. We obtain the quantum Langevin equations given by jl^j 

N S 

MQ k {t) = -DQ k (t)-K[2Q k (t)-Q k ^{t)-Q k+1 {t))-M^CkiQi(t) 

N s ft N s 

- E / ^ " *')&(*') dt ' - X>**(*)Q<(o) + a(*) + f(t), (9) 

e=i ^° i=\ 

with 

Mi kt = T(c kn 5 M -^), (10) 

Nb / ^ \ 

7w(*) = ( kn ~n J COSWnt, (11) 

Cfc(*) = y^c fcn I g n (0) cosu) n t + -4 — sinu> n £). (12) 

n=l V U « / 

Here £^ denotes the additional interaction between k and £th particles in the system induced 
by couplings {c kn }, lkn(t) stands for the memory kernel and Q k is the stochastic force. By 
using averages over initial values of q n (0) and q n (0), 

(mrtq n (0) 2 )B = m(qn(0f) B = coth , (13) 

we obtain the fluctuation-dissipation relation, 

l(C k (tMt') + Q(t')Ut)) B = E(^) (^)coth(^cos^(t-t'), (14) 

-)• fc B T 7w (t - for -+ 0, (15) 

where (•)# expresses the average over initial states of the bath. £ ke in Eq. (fTUj) denotes 
a shift of oscillator frequency due to an introduced coupling, and it vanishes if we adopt 



c n = mou^ f° r Ns = 1 [31]. In the case of Ns ^ 1, however, it is impossible to choose {ck n } 
such as £h = for all pairs of (k, £), then Qk is inevitably coupled with Qi (£ ^ k). Because 
of these couplings between HOs, the iVs-body system cannot be simply regarded as a sum 
of systems with Ng = 1. Although Eqs. f[9"j)- ([T2"j) are formally exact, it is difficult to solve 
iVs-coupled integrodifferential equations. 



B. The canonical transformation 



In order to obtain a tractable Langevin equation, we apply the canonical transformation 
to the model Hamiltonian. We assume that Ns is even without a loss of generality. Imposing 
a periodic boundary condition, 



QN s +k — Qk, PN s +k — Pk, 



(16) 



we employ the canonical transformation [20| . 

JV s /2-l 



J(2wks/N s ) 



Qs 



s=-N s /2 
N s /2-l 



p _ 1 M2-Kks/N s ) p 



(17) 
(18) 



-N s /2 



Note that the boundary condition is satisfied in Eqs. (|T7|) and ( f!8|) and that the set 
{(1/^Ns) e i( - 2nk / Ns ^ s } is orthogonal and complete in a periodic domain of the oscillator 
label k 2QJ]. By the canonical transformation, Hs in Eq. ((3]) becomes 

N s /2-l 



He 
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p;p s (D + Mn 2 s )Q* s Q s 

2M 2 



'N s Q f(t), 
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(20) 



where the commutation relations: 



iH5 ss >, 



[Q s ,Qs>] = [P s ,Ps>} = o, 



(21) 



hold with Q* = Q_ s and P* = P_ 
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For a simplicity of our calculation, we assume an identical frequency bath 

u n = u , c kn = c. (22) 

We furthermore assume that N B is even, imposing the periodic boundary condition given 
by 



QNb+u — ?n> PNb+u — Pn- 



We apply the canonical transformation 0, Q , 



N B /2-l 



1 °' 

7= y 



r=-N B /2 
Nb/2-1 
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; i(27rnr/Af s ) 



(23) 



(24) 



(25) 



r=-N B /2 



to the bath with the periodic condition given by Eq. ( 1231) . The bath Hamiltonian H B in 



Eqs. (jSJ) becomes 



25] 



M-/2-1 , „ 2 ~* ~ 



r=-N B /2 



2m 



where the commutation relations: 

[q r ,Pr>] = itfirr', [<?r, <?r'] = [Pr 3 2V] = 0, 



(26) 



(27) 



hold with g* = g_ r and p* = p- r . By canonical transformations given by Eqs. (fl71) . (fl8|) . 
d2H> and (ESD, tfj in Eq. © becomes 



s=-N s /2 r=-N B /2 



(28) 



Summing up Eqs. f lT9|) . fl26|) and (1281) . we obtain the total Hamiltonian expressed by 

H = H + H' S + H B , (29) 

where 
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with 

MQ 2 S = D + 4K sin 2 (j^-j + cN B for s = -fs ■ • ■, ^ - 1, (33) 
mCj 2 r = mwl + cN s for r = -^f , • • ^ - 1. (34) 

It is noted that H expresses the Hamiltonian for a uniform mode with s = u = and that 
a summation over s (r) in the H' s (H' B ) is excluded for s = (r = 0). 



C. Eigenfrequencies with f(t) =0 

Eigenfrequencies of the system-plus-bath with f(t) =0 may be obtained when we diag- 
onalize Hq given by Eq. ( )30|) . We employ the canonical transformation given by 

Q = M- 1 / 2 ^ cos 9 + X 2 sin9), P = M 1/2 (Y 1 cos 9 + Y 2 sin 9), (35) 
q = m _1/2 (-Xi sin^ + X 2 cos 9), p = m l/2 {-Yx sin9 + Y 2 cosfl), (36) 

where Yj = Xj and their commutation relations are given by 



[X^] = iH5 ij7 [XitXj] = [Yi,Yj] = for i,j = 1,2. 
The canonical transformation yields the diagonalized Hamiltonian given by 



H = H + H s + H B , 



with 



Ho 
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2 
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where H' s 
4>\ and 02 
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(39) 
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and ifg are given by Eqs. (131]) and ( 132|) . respectively. With the use of Eq. p0|) . 
are alternatively expressed by 
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where + (— ) of a double sign is applied to ((pf)- 

In the equilibrium state with f(t) = 0, Eqs. fl3Tl) . fl32|) and fj43|) yield eigenfrequencies of 
{z/j} (z = 1 to N s + N B ) for H given by 
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In the limit of c = 0, eigenfrequencies become 
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given by Eq. ( 1431) correspond to v + and z/_, respective 



eigenfrequencies, the system energy Eg is given by 

d\nZ s 



14[ | obtained the same eigenfrequencies by an alternative method: <pi and <p2 

y, in Ref . [l4j . With the use of these 
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Tr and Tre denoting a full trace over all variables and a partial trace over bath variables, 
respectively 



D. Responses to external forces 

1. Driven quantum harmonic oscillators 



Quantum HOs driven by an external force have been discussed in Refs. 2jJ-|23). It has 
been shown that the average position of a quantum HO is expressed by an equation of 
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motion of relevant classical HO 2lN23| as follows. The Hamiltonian of a single HO with 
mass m and oscillating frequency cuq driven by a force F(t) is given by 2ll423| 



H 



i 11 
p muZx _,, , 

— + xF(t), 

2m 2 w 



(49) 



for which the Schrodinger equation is expressed by 

h 2 d 2 



2m dx' 



muj 2 x 2 . 
+ °- xF(t) 



<9$(x,t) 
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(50) 



By using a unitary transformation, we may obtain a solution of $(x,t) expressed by [24j 



$ n (x,t) = n (a; - w(t)) exp 



mw(i - w(t)) -E n t+ I L(t') dt' 

o 



(51) 



with 



L(t) = -mw 2 — -muj 2 w 2 + wF(t), 



Jtloq ( n + 



for n = 0,1,2 



(52) 
(53) 



Here 0„(x) and E n are wavefunction and eigenvalue, respectively, of the Schrodinger equation 
with F(t) = in Eq. ( H9l) . and io(t) obeys an equation of motion for a classical driven HO, 



mw(t) + mujQwft) = F(t). 



(54) 



Equation (IBTjl shows that the center of a wave packet moves with w(t). It implies that an 



average of time-dependent position is given by 



21 



x(t) = w(t), 



(55) 



where an overline denotes the quantum average and w(t) is a solution of Eq. (j54p . This is 
consistent with Ehrenfest's theorem. 



2. Open quantum system of harmonic oscillators 



In order to study the response of the open quantum HO under consideration, it is nec- 
essary to pursuit equations of classical motions after Husimi's method 
(I30p and ( 135]) . the total Hamiltonian with f(t)^0 becomes 



21 



23J. From Eqs. 



H = H + H' s + H' B , 



(56) 
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with 

H = H 01 + H 02 , (57) 
Hoi = T + ^-/W Xl/(t)c0 ^' (58) 



H02 = f + ^L- ) J^X 2 f(t) S mB, (59) 

where H' s and H' B are given by Eqs. ( 13~TT) and ( 1321) . respectively, and <p\ and 02 are given by 
Eqs. (141 p and (l4"2"j) . Hamiltonians H i and i?02 i n Eqs. ( 15 8 p and ( I59p express HOs driven 



by forces of yjN s /M f(t) cos6> and ^N s /M f(t) sin6>, respectively From H' s in Eq. (1ST]) , 
equations of motion for Q s with s ^ are given by 

MQ S = -MCl 2 s Q s fors^O, (60) 

while Eqs. (157|) and (158]) lead to those for s = 0, X\ and A 2 , given by 



X, -ofA 1 + \/^/(f)cos0, (61) 



X 2 = -<^X 2 + ^f[t)mne. (62) 

A solution for Qo(^) m ay be evaluated from solutions of X±(t) and A 2 (t) with the canonical 
transformation given by Eq. ( 1351) . 

After some manipulations, quantum-averaged solutions of Q s are given by 

Q s (t) = Q s (0) cos ft s t + sin Q s t for s ^ 0, (63) 

Q (t) = Q (0)A Q (t) + P (0)A P (t) + %(0)B q (t) + p (0)B p (t) + $(t) for s = 0, (64) 

with 

2 

= ^Oicos0it, (65) 
, / x 1 v— ^ a i sin 0jt 

= (66) 

i=i ri 

-B g (i) = ~yj^ cos 6 1 sin 6* (cos 0it — cos 02^), (67) 
-d-dU = ; cosflsmfe 1 , 68 

$ W = IT E (^) J^mMt-t')f(t')dt' } (69) 
ai = 1 — a 2 = cos 2 0, (70) 
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where Q s (0), P s (0), <? s (0) and p s (0) denote initial states. The response of the total output 
averaged over initial states is given by 

R(t) = (J2Q k (t)\ =V^(5o(*)) fl . (71) 
\ k I o 

E(|)/ Sin ^(*-*W)< (72) 



M 



where we employ the relations given by 

4(0)) q = (P.(P)) Q = (?.(0)> = (P.(0)> = 0, (73) 

the bracket (-) expressing an average over initial states. Equation (172)) leads to the suscep- 
tibility, 

*<«> = § ('4) 

i=i ^* 

whose Fourier transformation is given by 

* (w) = ^£m^)' (75) 

with poles at u — ±</>j. 

It should be noted that R(t) in Eq. (I72"|) is generally not proportional to ATg except for 
the c = case because <pi and depend on Ns as shown in Eqs. f HTj) . and (jTDjl . This 
point will be shortly demonstrated in numerical model calculations for sinusoidal and step 
forces in the following. 

A. Sinusoidal forces 

We apply a periodic monochromatic force, 

fit) = gsinut, (76) 

where u and g stand for the frequency and magnitude, respectively, of the force. Equations 
(E2D and (J76D yield 

du\ f N sg\ v^«i sin wt-w sin 
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FIG. 1: (Color online) Eigenfrequencies v% of HO systems with Ng = 10 subjected to a bath with 
N B = 100 for c = 0.0 (open circles) and c = 10.0 (filled circles) (D = K = M = m = 1.0 and 
= 1.0), solid and dashed curves being plotted only for a guide of the eye. 



In the resonant case of u ■ 
'N s g 



R(t) 



M 



i, R(t) is given by 

di (sin ut — ut cos ut) a 2 (02 sin ut — u sin cf> 2 t) 



2u 2 



UT 



(78) 



Expressions of R(t) in the resonance cases of w = 02 an d u = 0i = 02 are similarly given. 
In the limit of c = where 0i = Qq, 02 = ^o, = 0.0, a\ = 1.0 and a 2 = 0.0, Eq. ( 177j) 
reduces to 



M 



f2o si n cot — oj sin f2 t 



for c = and w/flc 



^(^-w 2 ) 

which expresses the response of a HO isolated from a bath. 

We have performed numerical model calculations, choosing a coupling fl3 |. 

c 



(79) 



(80) 



such that the interaction term in Eq. §6§ including summations over Ylk=i an d Enfi yield 
finite contributions even in the limits of Ng — > oo and/or Nb — > oo. We have adopted 
parameters of D — K — M — m — Uq — 1.0 for a given system-plus-bath. 

Figure [1] shows eigenfrequencies Ui for c = 0.0 (open circles) and c = 10.0 (filled circles) 
of a HO system (Ng = 10) subjected to a bath (iV s = 100). Eigenfrequencies z/j for 
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FIG. 2: (Color online) Responses of R(t)/Ns of HO systems with N$ = 1 (chain curves), 2 
(dashed curves), 10 (dotted curves) and 20 (solid curves) for (a) Co = 1.0 and (b) Co = 10.0 to 
an applied sinusoidal force with u = 0.5 and g = 1.0 (D = K = M = m = 1.0, ojq = 1.0 and 
N B = 100). 

1 < i < 10 show a dispersion relation of the HO system while those for 11 < < 110 of 
the bath are almost constant. For cq = 0.0, we obtain = 1.0 and ujq = 1.0. When the 
system-bath coupling of Co = 10.0 is introduced, they become 1.414 and 1.048, respectively, 
which lead to <pi = 1.449 and <p2 = 1-0. 

Figure a) shows responses of R(t)/N s to a sinusoidal force with uj = 0.5 and g = 1.0 
of HO systems with N$ = 1, 2, 10 and 20 coupled to Nb = 100 baths with a coupling 
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<l>2 (C =l)=4>2 (CQ=10) 




FIG. 3: (Color online) N$ dependences of (a) fa, and (b) 8 and aj (i = 1,2) for cq = 1.0 (dashed 
curves) and 10.0 (solid curves) (D = K = M = m = ujQ = 1.0 and Nb = 100). 6 and ai in (b) are 
independent of Co for the adopted parameters (see the text). 

of Co = 1.0. Results of R(t)/Ns are almost the same independently of N$, although some 
discrepancies among the four results are realized at t > 50. These discrepancies become 
more evident for a larger coupling of c = 10.0, whose results are shown in Fig. [2(b). These 
results clearly suggest 

R(t; N s ) = N s R(t; 1) for c = 0, (81) 
^ N s R(t; 1) for c ^ 0. (82) 
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In order to elucidate Ns and Co dependences of R(t)/Ns, we show in Figs. [3j 0i, 6 and 
Oj (i — 1, 2) as a function of iVs for Co = 1.0 (dashed curves) and 10.0 (solid curves). Figure 
El^a) shows that with increasing Ns, 0i is slightly increased while 02 is constant. We note 
in Fig. E(b) that an increase of Ng yields an increase in 9, by which 02 is increased but 
ai is decreased. For adopted parameters, 9, ai and a 2 are independent of c because the 
denominator of Eq. (H0|) becomes — Uq = (Nb — Ns)c whose c is cancelled out by that 
in its numerator. With increasing Ns, a contribution from a lower eigenfrequency of 02 is 
increased. The effect of the system-bath coupling for cq = 10.0 is more significant than that 
for Co = 1.0 because the difference of 0i — 02 in the former is larger than that in the latter: 
if 0i = 02 results are independent of Oj (and then N s ). 

B. Step forces 

Next we apply a step force given by 



where Q(x) stands for the Heaviside function and t s is the starting time of a force with a 
magnitude of g. The averaged output is given by 



FigureH^a) shows R(t)/Ns for a step force with t s = 10.0 and g = 1.0 of HO systems with 
N s = 1, 2, 10 and 20 coupled to N B = 100 baths with a coupling of c = 1.0 (D = 1.0 and 
Uo = 1.0). Result of R(t)/Ns for Ns > 2 are almost the same as that for Ns = 1. However, 
when the interaction is increased to cq = 10.0, the discrepancy between results of Ns = 1 
and Ns > 2 become evident. Fig. H^b) shows similar plots but with stronger coupling of 
Co = 10.0, for which shape and magnitude of R(t)/Ns are significantly modified for Ns > 2. 

III. DISCUSSION 

A. Responses to space- and time-dependent forces 

It is interesting to calculate responses to a space- and time- dependent force which yields 



f(t) = g e(t s -t) 



(83) 




(84) 



H f in Eq. © 



H f = -f(t)S(u), 



(85) 
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FIG. 4: (Color online) Responses of R(t)/Ng of HO systems with N$ = 1 (chain curves), 2 
(dashed curves), 10 (dotted curves) and 20 (solid curves) for (a) cq = 1.0 and (b) Co = 10.0 to an 
applied step force with t s = 10.0 and g = 1.0 (B = K = M = m = 1.0, oj = 1.0 and N B = 100). 
Results for all N$ in (a) are indistinguishable. 

with 

N B 

S( u ) = J^Q k e l27Tku/Ns for u e {-fs + 1, ^ - 1}. (86) 

k=l 



16 
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FIG. 5: (Color online) Responses of R[t)/Ns of an isolated HO system (co = 0.0) to an applied 
sinusoidal force with u = 0.5 and g = 1.0 for various u (Ns = 12, Nb = 100 and K = D = M = 
m = ojq = 1.0). Results for u = —6.0, —4.0 and —2.0 are shifted by 6.0, 4.0 and 2.0, respectively, 
for clarity of the figures. 



Here the wavevector u is, for example, u = and u 
forces, respectively, for which S(u) is represented by 

N S 

fe=i 

N S 



-Ng/2 for uniform and staggered 



S(u) 



for u = 0, 



fc=i 



k e 



-iirk 



for m 



TVs. 
2 ' 



(87) 



(88) 



The mode with m ^ does not couple with s = mode which couples with bath as 
mentioned in the preceding subsection II D. Equations of motion for Q s with s ^ are 
independent of degrees of freedom in a bath and they are given by 

MQ U = -Mtt 2 u Q u + f(t) for s ^ and s = u ^ 0, (89) 
MQ S = -Mtt 2 s Q s for s ^ and s ^ u ^ 0. (90) 

The response to applied force with u 0) is given by 

'N s g\ f f sin n u (t-t')f(t') 



R(t) 



M 



(91) 
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FIG. 6: (Color online) Responses of R(t)/Ns of HO systems coupled with Co = 0.0 (dashed curves), 
Co = 5.0 (chain curves) and Co = 10.0 (solid curves) to an applied sinusoidal force with ui = 0.5 and 
g = 1.0 for (a) u = 0.0 and (b) u = -6.0 (N s = 12, Nb = 100 and K = B = M = m = oJo = 1.0). 
Results for co = 0.0, 5.0 in (a) are shifted by 4.0 and 2.0, respectively, and those for co = 0.0, 5.0 
in (b) are similarly shifted by 2.0 and 1.0, for clarity of the figures. 



(92) 



which becomes for sinusoidal force [Eq. ( 1TB]) ]. 

R(t) - ( Nsg \ ( ^ u s^ut - usinttj 

u " V m ) \ n u (ni-^) 

In the limit of c = 0.0, R(t) is given by Eqs. <^2) and (JH]) with tt u = O u . The effect of 
finite coupling is realized by a change in Q u as given by Eq. (I33I) . Note that the response 
to applied force with u = has been studied in subsection II D [Eq. f!72]) ]. 



18 



We present model calculations for sinusoidal forces with u = 0.5 and g = 1.0 in Eq. (!76|) 
for N s = 12, N B = 100, K = B = M = m = ui = 1.0. Figure E shows R(t)/N s for 
isolated systems (c = 0.0) with u = 0.0, -2.0, -4.0 and -6.0. Magnitudes of R(t)/N s 
become smaller for larger \u\. Figure a) and EJb) show R(t)/N$ for uniform (u = 0.0) 
and staggered forces (u = —6.0), respectively, with couplings of c = 0.0 (dashed curve), 5.0 
(chain curve) and 10.0 (solid curve). Comparing Fig. Mb) with Fig. E^a), we notice that 
an effect of couplings for staggered forces is less effective than that for uniform forces. 



B. Spatial correlation 



Employing eigenfrequencies for f(t) = obtained in subsection II C, we may calculate 
the spatial correlation between and Qk+ m , 



N S 



Tm = ^ (QkQk+m) 



k=l 

Ns/2-1 



\Q:Qs)e- i27Tms/Ns , 

s=-N s /2 



(93) 
(94) 



with (Q* s Qs) evaluated by [19] 



QIQ? 



dlnZc 



(95) 



,(3MQ S J dQ s 

where the bracket () denotes the average over H and Z$ is given by Eq. fj46|) . T m with 
m = expresses a (summed) variance of Qk'- r = ^2k=i(Qt) ■ After some manipulations 
with the use of the diagonalized Hamiltonian given by Eq. (138 p . we obtain 



q* s Qs 



h ( (3hn s 

-=— coth 



2Ma 



V 2 



-^Vcoth^ 



with 



<90i 


fio 


<9fi 


" 20! 


(902 


fi 


dfl 


"~ 20 2 



^^0^ 90i 
<9fir 



O 2 , ! 2 



for s ^ 0, 



for s = 0, 



(96) 
(97) 



(fig _ w 2 ) 2 + 4]\[ s N B c 2 / Mm 

6 2 , , 2 
"o ~ 

(fig - c^) 2 + 4iV s iV B c 2 /Mm 



(98) 



(99) 
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where 4>i an d 02 are given by Eqs. (|40j) and (HIT) . Substituting Eqs. (1961) and (1971) into Eq. 
(??), we obtain r m , 



+ 



Ar s /2-l 

E 

s=-N s /2 

h 

2MQn 



— coth — — 

\ 2 



coth 




—i2irms/N s 



For T 



and T — )• oo, T m becomes 



2Ma 



jV s /2-l 

E 

s=-7V s /2 
Ms/2-1 / , _ x 



^—i2irms/N s _|_ 



■i2irms/N B 



+ 



2Mfi 



i=l 



coth 



90i_ 



o 



E 

4 = 1 



9 In c; 



1 



(100) 



for T = 0, (101) 



for T ^ oo. (102) 



In the case of uncoupled, isolated system with K = 0.0 and c = 0.0, T m is given by 



5 



m0 



2Mn f 



coth 



for if = 0.0 and c = 0.0, 



(103) 



which is proportional to Ng and which vanishes for m > 1. It is, however, not the case for 
if 0.0 or 0.0. Indeed in the case of K ^ 0.0, T m is finite for m > 1 because of direct 
particle-particle couplings of K and indirect couplings of —CkiCn n jmuj\ in the second term of 
Eq. ( flOl) . Even when if = 0.0, T m with c 7^ 0.0 remains finite with a small negative value. 

Figure [7] shows the temperature dependence of T m (T)/N s for m = (solid curve), 1 
(dashed curve) and 2 (chain curve) of HO systems with N$ = 10, Nb = 100, K = D = M = 
m = ujq = 1.0 and c = 10.0. T m (T) is finite at T = 0, and at T — > 00 it is proportional to 
temperature, as Eqs. (11011) and ( 11021) show. Magnitude of T m is smaller for a larger m. The 
dotted curve expresses Co (= (QqQo)) which is larger than r because flo < Q s with s 7^ 0. 

iVs dependences of T Q (T)/Ns at ksT/hwo = 0.0 and 10.0 are shown in Fig. El^a) and 
Mjo), respectively, for c = 0.0 (open circles), 1.0 (filled square) and 10.0 (filled circles). For 
c = 0.0, r (0) is proportional to N s as expected. However, when the system-bath coupling 
is introduced, T is not proportional to Ng as shown in Fig. |HJ This is realized not only at 
zero temperature but also at high temperature. 

Even when external forces are applied, the spatial correlation is not modified, which is 
the characteristics of the open system with the linear system-bath coupling. In the open 
system with the nonlinear system-bath coupling, the spatial correlation is modified by an 



applied force 



26]. 
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FIG. 7: (Color online) The temperature dependence of T m (T)/Ns for m = (solid curve), m = 1 
(dashed curve), m = 2 (chain curve) and Co (= {QqQo)) (dotted curve) for a HO system (N$ = 10, 
iVB = 100, K = D = M = m = ojQ = 1.0 and cq = 10.0) 

IV. CONCLUSION 



Responses of open small quantum systems described by the (N s + N B ) model [13|, [14 1 



have been studied. By using double canonical transformations mentioned in subsections II B 
and II C, we obtain the diagonalized Hamiltonian, from which the response to applied forces 
is obtained with the use of Husimi's method for a driven quantum HO [21]. The response 
to a uniform force given by Eq. (jlj) is generally not proportional to N$ against our implicit 
expectation. This nonlinear response is consistent with the system specific heat in open 
small quantum systems previously discussed in Ref. and it is realized also in spatial 

correlation T m not only at low temperatures but also at high temperatures. These facts 
show an importance of taking account of finite N$ in discussing open quantum and classical 
systems. It would be interesting to examine the obtained non-linearly by experiments for 
open small systems. 



21 




FIG. 8: (Color online) The Ns dependence of To(T)/Ns at (a) ksT/foLOQ = 0.0 and (b) hsT/tno® = 
10.0 for Co = 0.0 (dashed curve), 1.0 (chain curve) and 10.0 (solid curve) of a HO system {Ng = 10, 
N B = 100, K = D = M = m = uj = 1.0). 
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